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ABSTRACT 

Connecting cosmological simulations to real-world observational programs is often complicated by 
a mismatch in geometry: while surveys often cover highly irregular cosmological volumes, simulations 
are customarily performed in a periodic cube. We describe a technique to remap this cube into 
elongated box-like shapes that are more useful for many applications. The remappings are one-to-one, 
volume-preserving, keep local structures intact, and involve minimal computational overhead. 
Subject headings: methods: iV-body simulations — cosmology: large-scale structure of universe 



1. INTRODUCTION 

Numerical simulations have become an indispensable 
tool in modern cosmological research, used for investigat- 
ing the interplay of complex physical processes, studying 
regimes of a theory which cannot be attacked analyti- 
cally, generating high precision predictions for cosmolog- 
ical models, and making mock catalogs for the interpre- 
tation and analysis of observations. Such simulations tra- 
ditionally evolve the matter distribution in periodic cu- 
bical volumes, which neatly allow them to approach the 
homogeneous Friedmann solution on large scales. The 
use of a periodic volume also allows the long-range force 
to be easily computed by fast Fourier transform methods 
in many popular algorithms (e.g. particle-mesh, particle- 
particle-particle-mesh or tree-particle-mesh algorithms). 

Surveys, on the other hand, often cover cosmological 
volumes that are far from cubical in shape, and making 
a mock catalog which includes the full geometrical con- 
straints of the observations is difficult. One approach is 
to simulate a sufficiently large volume that the survey 
can be embedded directly within the cube, but this often 
means large parts of the computational domain are un- 
used. An alternative is to trace through the cube across 
periodic boundaries so as to generate the desired depth, 
with various rules for avoiding replication or double- 
counting of the volume (if desired). A third approach 
is to run a simulation in a non-cubical geometry. If the 
side lengths are highly disproportionate this can lead to 
its own numerical issues, and in addition it makes it dif- 
ficult to reuse a given simulation for many applications. 

In this paper we present a new solution to this prob- 
lem which allows one to embed a hypothetical survey 
volume inside a cosmological simulation while limiting 
wasted volume and artificial correlations. The method 
is based on the simple observations that (1) a cube with 
periodic boundary conditions is equivalent to an infinite 
3-dimensional space with discrete translational symme- 
try, and (2) the primitive cell for such a space need not 
be a cube. We show that one may take the primitive 
cell to be a cuboicQ of dimensions Li x L2 x L3 for a 
discrete but large choice of values Li > L2 > L3. The 

^ To be precise, by cuboid we mean here a rectangular paral- 
lelepiped, i.e. a parallelepiped whose six faces are all rectangles 
meeting at right angles. For practical purposes, a cuboid is simply 
a box that is not a cube. 
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Figure 1. Possible dimensions Li X L2 X L3 for cuboid remap- 
pings, subject to the conditions Li > L2 > L3, L1L2L3 = 1, and 
Li < 7. The choice used for the Stripe 82 remapping described in 
Section |4] is indicated by a red circle. 

possible choices, subject to the constraint L1L2L3 — 1, 
are illustrated in Figure [1] 

Our approach leads to a one-to-one remapping of the 
periodic cube which keeps structures intact, does not 
map originally distant pieces of the survey close to- 
gether, and uses no piece of the volume more than 
once. It complements existing techniques fo r gen- 
erating mock observation s (e.g. iBlaizot et al.l 120051 : 
iKitzbichler fc White! I2007|) or for ray-tracing through 



simulations (e .g. 'White fc H d[2000|: IVale fc White! [2001 
iHilbert et all 12009: FuUana etal.1 12010[ ). Though not 



ideal in all cases, our remapping procedure may be seen 
as a general-purpose alternative that neatly skirts many 
of the complications involved with previous methods. 

We begin in Section [2] with a mathematical description 
of the remapping, and explain what choices of dimensions 
are possible. In Section [3] we describe how to implement 
this remapping numerically. We present a few useful ex- 
amples in Section 21 and conclude in Section [5] with a 
discussion of some of the advantages and limitations of 
this method. 

2. MATHEMATICAL DESCRIPTION 
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Consider a unit cube with periodic boundary condi- 
tions. By tihng copies of this cube in all directions, each 
point (x, y, z) G corresponds to a point in the canon- 
ical unit cube [0, 1]^ via the map 

(x, y, z) H> {x mod 1, y mod 1, z mod 1), (1) 

where "a; mod 1" is just the fractional part of a; if a; > 0, 
or one minus this fractional part for x < 0. Since we can 
identify each point in M.^ with a point in [0, 1]'^, any re- 
gion in space corresponds to a sampling of the unit cube. 
We are primarily interested in bijective samplings, i.e. re- 
gions that according to the above equivalence cover each 
point of the unit cube once and only once. We call such 
a sampling a remapping of the unit cube. Equivalently, 
a remapping may be thought of as a partitioning of the 
unit cube into disjoint regions, which are then translated 
by integer offsets and glued back together along periodic 
boundaries. 

A general class of such remappings may be constructed 
by applying shear transformations to the unit cube. We 
first illustrate the idea in 2D, where the transformations 
may be easily visualized, and then state the appropriate 
generalization to 3D. Start with a continuous field laid 
down within the unit square, which may be thought of as 
a parallelogram spanned by the vectors Ui = (1,0) and 
U2 — (0,1). We suppose the field to be periodic, and 
tile copies throughout the plane [see Figure [2ja)]. Now 
imagine taking the unit square and shearing it along its 
top edge, leaving the underlying field in place [Figure 
[DJb)]. The result is a new parallelogram defined by the 
vectors 

Ml = (1,0), M2 = M2 + mMi = (to, 1), (2) 

where m is some real scalar that controls the extent of the 
shear. Since shear transformations are area-preserving, 
this parallelogram has unit area. In fact, since Mi is 
a lattice vector (i.e. a vector of translational invariance), 
this parallelogram covers each point of our field once and 
only once, and hence defines a valid remapping. 

Now observe that if m is an integer, then Mj — (to, 1) 
will also be a lattice vector. In that case we again have 
a parallelogram with edge vectors Ui and u'2 that are 
both lattice vectors. We can now shear this parallel- 
ogram along its right edge [Figure [SJc)], giving a new 
parallelogram with edge vectors 

m'i = Ml + nu'2 = (1 + mn, n), u'2 — (to, 1). (3) 

Once again this parallelogram covers each point of our 
field exactly once, and if n is an integer then u[ and u'2 
are both lattice vectors. We may repeat this process, 
applying integer shears alternately to the top and right 
edges of the parallelogram. In general, if Mi and M2 are 
any integer-valued 2D vectors that span a parallelogram 
of unit area, they can be obtained by such a sequence of 
shear transformations. This condition is simply 

U11U22 - W12W21 = det ( ) = 1, (4) 

\U21 U22 J 

i.e. Ml and M2 are the rows of an invertible integer-valued 
matrix. 

While remapping the unit square into a parallelogram 
may be useful in certain cases, parallelograms are gen- 
erally too awkward to be useful. Instead, given lattice 



vectors Mi and M2 spanning unit area, we may apply one 
last shear to "square up" the parallelogram into a rect- 
angle [Figure [2](d)]. Explicitly, we let 

ei = Ml, e2=M2 + Q;Mi, (5) 

and choose a so that ei and 62 are orthogonal. The 
rectangle defined by ei and 62 now covers each point of 
the unit square exactly once. 

The generalization to 3D is straightforward. By apply- 
ing integer shears to the faces of the cube, we may obtain 
any parallelepiped with edges given by integer vectors Mi , 
M2, M3 satisfying 

/Uii U12 Uiz\ 

det U21 U22 U23 =1- (6) 

\M31 W32 U33/ 

Again this space of possibilities corresponds to the space 
of invertible, integer-valued 3x3 matrices. We apply 
two final shears to square up this parallelepiped into a 
cuboid, by choosing coefficients a, /?, and 7 such that 

ei = Ml, (7) 
62 = M2 + aMi, (8) 
ea = M3 + /?Mi + 7M2, (9) 

are mutually orthogonal. This gives a remapping of the 
unit cube into a cuboid with side lengths Li = |ei|, sat- 
isfying L1L2L3 = 1. Moreover since ei is still a lattice 
vector, this cuboid is periodic across the face perpendic- 
ular to this vector. 

Invertible, integer-valued 3x3 matrices can be gener- 
ated easily by brute force computation, and each such 
matrix leads to a cuboid remapping of the unit cube. 
The allowable dimensions Li x L2 x L3 for these cuboids 
are illustrated in Figure [U where L3 is given implicitly 
by the unit volume condition L1L2L3 = 1. 

3. NUMERICAL ALGORITHM 

The goal of our remapping procedure is to provide an 
explicit bijective map between the unit cube and a cuboid 
of dimensions Li x L2 x L3. We will refer to points in 
the imit cube by their simulation coordinates x € [0, 1]"^, 
and their remapped positions in a canonical, axis-aligned 
cuboid by remapped coordinates r G [0, Li] x [0, L2] x 

[0,^3]. 

The reverse map r ^ x is simple. The edge vectors 
ei, 62, 63 discussed previously describe how to embed an 
oriented cuboid within an infinite tiling of the unit cube. 
Let fii = Bi/Li be unit vectors along these edges. Then 
given (ri,r2,r3) e [0,Li] x [0,^2] x [0,^3], the point 

p = nhi + r2fi2 + r^ns (10) 

lies within this oriented cuboid, and this maps to a point 
in the unit cube according to Eq. ([T]). 

The forward map a; 1-^ r is slightly more complicated. 
Consider again the oriented cuboid with edge vectors 
61,62,63 embedded within an infinite tiling of the unit 
cube. Each tile (i.e. each replication of the unit cube) 
may be labeled naturally by an integer triplet m € li^, 
where the canonical unit cube has m = (0,0,0). The 
intersection of the cuboid with a tile is called a cell, of 
which only a finite number will be non-empty. Each non- 
empty cell is uniquely labeled by the triplet m. [The four 
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Figure 2. Panels (a)-(d) show a sequence of shear transformations that results in a remapping of the unit square into a rectangle of 
dimensions \f% X l/\/5, as described in the text. Panel (e) highlights the edge vectors ei, 62, and the four cells that comprise this rectangle. 
Note that each point of the original unit square is covered exactly once in the final remapping. 

cells of the 2D example from the previous section are in- 
dicated in Figure [Sfe) .] 

As a geometrical object, each cell is just a convex poly- 
hedron bounded by 12 planes: the 6 faces of the tile and 
the 6 faces of the cuboid. A plane may be parametrized 
by real numbers (a, &, c, d), so that a point {x,y,z) lies 
inside, outside, or on the plane depending on whether 
the quantity ax + by + cz + d is less than, greater than, or 
equal to zero. Thus to test whether a point belongs to a 
cell, we need only check if it lies inside all the planes that 
bound it. When translated spatially by a displacement 
— m, each cell defines a region within the unit cube. The 
collection of all such cells defines a partitioning of the 
unit cube, with each point in [0, 1]'^ being covered by ex- 
actly one cell. Our algorithm for the forward map x t-^ r 
then may be summarized as follows: 



1. determine which cell contains the point x (using 
point-plane tests) 

2. let p = a; -I- m be the corresponding point in the 
oriented cuboid 

3. define r e [0, Li] x [0, L2] x [0, L3] by ri = p ■ hi 

4. EXAMPLES AND APPLICATIONS 

The mappings described above can be used in myriad 
ways. The algorithm is fast enough that it can be used for 
on-the-fly analysis while a simulation is running, it intro- 
duces little overhead in walking a merger tree (allowing 
complex light-cone outputs to be built during e.g. the 
running of a semi- analytic galaxy formation code), or it 
can be applied in post-processing to a static time output 
of a simulation to alter the geometry. In this section we 
give a (very) few illustrative examples. 

As it is so widely known we sh all use as our fiducia l 
volume the Millennium simulation (jSpringel et al.|[2i305D . 
which was performed in a cubical volume of side length 
500/i~^Mpc. We begin by asking how we could embed 
a very wide-angle survey, such as the equatorial stripe 
("Stripe 82") of the Sloan Digital Sky Surve}^!, in such 
a volume. Stripe 82 is 100° wide and 2.5° in height. 
Using the Millennium simulation cosmology, the volume 
within Stripe 82 out to z ~ 0.45 equals the total vol- 
ume within the Millennium simulation. At this depth 
the stripe is 1.2/i~^Gpc in the line-of-sight direction, 
1.9/i~^Gpc transverse but only 50/i~^Mpc deep. 

We can map the cubical volume of the simulation into 
this geometry using the transformation circled in red in 
Fig.[TJ A view of this transformation is given in Figure [3l 
where the regions within the cube are shown "unfolded" 

^ |http: //www. sdss . org| 




Figure 3. The regions in the unit cube, and how they "unfold" 
to produce the wide, thin geometry shown in Figure [4] 

to produce a long and wide, but thin, domain into which 
we can embed the Stripe 82 geometry. Figure 0] shows 
a light-cone produced from the outputs of the Millen- 
nium simulation, in the Stripe 82 geometry. Specifically 
we show all galaxies brighter than R = —20.5 (to avoid 
saturatin g the figure) and with < z < 0.45 from the 
catalog of iDe Lucia fc BlaizotI (|2007[ ). 

Another frequently encountered situation is a survey 
which is much longer in the line-of-sight direction than 
either of the (approximately equal) transverse directions. 
Let us consider for example a survey which aims to reach 
z ~ 3, or about 5 h~^Gpc in the Millennium simula- 
tion cosmology. Among the many possible choices at 
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Figure 4. A mock catalog from the Millennium simulation remapped into the equatorial stripe of the Sloan Digital Sky Survey. The 
"stripe" is a thin rectangle on the sky of angular dimensions 100° X 2.5°, shown here "from above" extending out to redshift z ~ 0.45, or 
r ~ 1.2/i~^Gpc comoving. We show only galaxies brighter than R = —20.5, to avoid saturating the figure. 



our disposal we find a remapping with {L\^Li^L-i) = 
(5025, 179, 139) /i~^Mpc, which could encompass a sur- 
vey of angular dimensions 2.0° x 1.6° out to 2 ~ 3.5. 
To probe even earlier epochs we could choose a remap- 
ping with (Li, L2, is) = (7106, 145, 121) /I'^Mpc, which 
allows a 1.2° X 1.0° survey out to z ~ 10. For these types 
of highly elongated geometries, the possible remappings 
are dense enough that a suitable choice may be found for 
almost any survey. 



5. DISCUSSION 

As surveys become increasingly complex and pow- 
erful and the questions we ask of them become in- 
creasingly sophisticated, mock catalogs and simula- 
tions which can mimic as much as possible the obser- 
vational no n-idealities become increasingly important. 
lAngulo fc~ Whitc (2009) have shown that simulations of 
one cosmology can be rescaled to approximate those of 
a different cosmology. We have introduced a remapping 
of periodic simulation cubes which allows one simulation 
to take on the characteristics of many different observa- 
tional geometries. The use of such techniques enhances 
the usefulness of cosmological simulations, which often 
involve a large investment of community resources. 

The methods we introduced in this paper lead to one- 
to-one remappings of the periodic cube which keep struc- 
tures intact, do not map originally distant pieces of the 
survey close together, and use no piece of the volume 
more than once. The remapping can be done extremely 
quickly, meaning it can be included in almost any anal- 
ysis tool with negligible overhead. 

Remapping the computational geometry is, however, 
not without its limitations. First and foremost, although 
the target geometry may have sides much longer than 
the original simulation, the structures will not contain 
the correct large scale power since it was missing from 
the simulation to begin with. This problem becomes less 



acute as the original simulation volume becomes a fairer 
representation of the Universe. Secondly, if the target ge- 
ometry is too thin it is possible for points which are far 
apart in the survey to come from points close together in 
the simulation volume, leading to spurious correlations. 
These are analogous to the artificial correlations between 
survey "sides" that occur when it is embedded in a peri- 
odic cube. Simply excluding a boundary layer from the 
remapped volume can tame such correlations. 

In addition to the description of the algorithm in 
this paper, we have made Python and C-l — h im- 
plementations of the remappings, along with fur- 
ther examples and animations, publicly available at 
,http : / /mwhite . berkeley . edu/BoxRemap , 
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